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Abstract 

The determination of the smuon mass at the Clic multi-TeV e + e~ linear col- 
lider has been studied for two CMSSM benchmarks. Results are given for both 
the analysis of the muon energy spectrum and the threshold scan method. The 
effects of detector resolution, beam-beam interactions and accelerator-induced 
backgrounds are discussed. The energy spectrum technique is also applied to 
the t — > tg process to determine the scalar top mass, in scenarios with the gluino 
lighter than the squarks. 



1 Introduction 



A multi-TeV e + e~ linear collider (LC) is expected to complete the investigation of 
the SUSY particle spectrum, in the case Supersymmetry is realised in nature and the 
particle partners exhibit a mass spectrum extending up to 1 TeV, and beyond. In this 
case, precise determinations of the masses and widths of these particles will represent 
a central part of the collider physics programme. Sets of SUSY benchmarks have 
recently been proposed, which take into account constraints from direct searches at 
Lep-2 and other colliders and data from b — > S7 and cosmology [2, 3]. Most of 
these benchmarks have the masses of the heaviest sparticles beyond the reach of the 
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Lhc and of a TeV-class e + e~ linear collider. The Lhc is expected to reveal the 
existence of Supersymmetry and perform a study of the sparticle properties in the 
vast majority of conceivable scenarios and a 0.5-1.0 TeV e + e~ collider can improve 
the accuracy in the study of the lightest SUSY states. However, higher energies are 
likely to be needed to access all the sparticles species and acquire the data necessary 
to clarify the nature of Supersymmetry breaking. 

A first mapping of the potential of a 3-5 TeV e + e~ collider, such as Clic [4], in 
comparison to the Lhc and to a lower energy e + e~ collider was performed, based on 
a simple assumption on the minimum observable product of production cross sec- 
tion and decay branching fraction into well reconstructable final states [2]. Beyond 
detection, it is important to assess the ability of Clic to perform accurate deter- 
minations of sparticle masses, widths and decay branching fractions, taking into 
account the details of beam-beam interactions and the resulting luminosity spec- 
trum and accelerator-induced backgrounds [1]. We study here the expected Clic 
accuracy on sfermion masses, using a realistic simulation. Results are presented for 
the Supersymmetric partner of the muon, fi^ and for the scalar top, £2- Section 2 
presents the benchmarks considered and the Clic parameters. There are two differ- 
ent techniques for the determination of the smuon mass: the analysis of the energy 
distribution of the resulting muon in the two body /II decay and the energy scan 
close the the production threshold. These are discussed in 3.1 and 3.2 respectively. 
In Section 4 we discuss the perspectives for extending the energy spectrum analysis 
to ti decays and Section 5 has the conclusions. 



2 SUSY Benchmark and Simulation 



The post-LEP CMSSM benchmark points E and H, proposed in [2], have been chosen 
for this study. Point E 1 is representative of the focus region, characterised by heavy 
sleptons and squarks. The fiL + LiL~ threshold is just accessible at y/s =3 TeV. The 
gluino (Mg=697 GeV) being lighter than the squarks, decays such as q — > qg are open 
for these parameters. It is interesting to observe that such inverted scenario, with the 
gluino lighter than the squarks, is also common to benchmarks inspired to effective 
field theories derived from strings [5]. The mass of the LSP is 120 GeV and those of 
the heavier neutralinos are in the range 200-300 GeV. Point H 2 is located at the end- 

1 m 1/2 = 300, m = 1500, tan/3 = 10, sign(n) = +, A = 
2 m 1/2 = 1500, m = 419, tan/3 = 20, signal) = +, A = 
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Table 1: Main properties of the sparticles considered in this analysis 



Point 




Mass 
(GeV) 


Width 

(GeV) 


a(e + e — > fiji) (fb) 
(fb) 


BR(/i - M X?) 


E 


ft-R 


1433 


21. 


0.03 


0.35 


E 




1427 


7.9 


0.02 


0.95 


H 




710 


0.06 


0.92 


1.00 


H 


fJ-L 


1150 


0.6 


0.59 


1.00 










a(e + e~ -► tt) (fb) 


BR(t -► tflr) 


E 


h 


994 


42. 


0.65 


0.20 


E 


*2 


1303 


76. 


0.12 


0.48 



point of the coannihilation tail at large and moderate tan f3. The heavy sparticle 
spectrum makes this a highly problematic point at the Lhc. However, at ^/s=3 TeV, 
smuons are accessible in e + e~ collisions and their masses can be determined through 
the muon energy distribution. The LSP mass, M o, is 665 GeV. Events have been 
simulated using the Pythia 6.2 Monte Carlo generator. Parameters have been 
tuned to reproduce sparticles masses and widths to a good accuracy. Initial state 
radiation (ISR) has been included. Muon energy distributions have been obtained 
using full Geant-3 simulation of a discrete silicon tracker consisting of twelve layers 
located from 3.5 cm up to 180 cm in radius and a solenoidal magnetic field B of 
either 4 or 6 T. Backgrounds and t decays have been processed using the Simdet 
parametric simulation [6], describing the detector response determined using Geant. 
The Clic luminosity spectrum has been folded to the production cross sections using 
the Calypso interface. 77 — > hadrons and parallel muon backgrounds have been 
overlayed to the generated signal and background events for the smuon study. 



3 Muon Mass Determination 



The study has been performed for e + e~ — ► — * A t+ XiA* - Xi- The three main 
sources of background, also leading to two muons plus missing energy, are i) e + e~ — > 
W + W~ — > n + ^Tv^u^, ii) e + e~ — > W + W~vv — > ^ a~v il y^v e v e and iii) e + e~ — > 
X1X2, X2X2 — /W + /U _ ^^XoXo- These backgrounds can be suppressed by requiring 
central production and decay kinematics compatible with those characteristic of 
smuon pair production. A multi-dimensional discriminant based on M^, M recoi i, 
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E 



missing 



, /ifj, Acolinearity, | cos ^Thrust I , E t and Ef iem has been applied. The signal 



efficiency is flat with the muon energy. 



3.1 The Energy Distribution method 



If the centre-of-mass energy, y/s, is significantly larger than twice the sparticle mass, 
Mn, this can be determined by an analysis of the energy spectrum of the muon, 
emitted in the two-body jl — ► XjjU decay (see Figure 3.1). The two end-points, E min 
and E max , of the spectrum are related to the jl and Xi masses and to the jl boost 
by: 



E, 



max I min 



2 v 



Xi 

Ml 



x (1 ± 



\ 



E 2 

beam 



from which either the smuon mass can be extracted, if M^o is already known, 



or both masses can be simultaneously fitted. This technique, already considered 




700 



680 - 



640 




100 200 300 400 500 600 700 800 900 1000 

P u (GeV) 



1110 1120 1130 1140 1150 1160 1170 1180 1190 
Smuon Mass (GeV) 



Figure 1: Left: Muon energy spectra in the decay jl — > fiXi f° r = 1150 GeV and 
M x o = 660 GeV obtained at \fs = 3 TeV, assuming the baseline Clic luminosity 
spectrum. Right: Accuracy on the determination of the mu and x± masses by a 
two-parameter fit to the muon energy distribution. The lines give the contours at 
1 a, 68% and 95% C.L. for 1 ab~ l of data at ^=3 TeV. 
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Table 2: Results of the 1-parameter x 2 &t to the muon energy distribution, obtained 
for different assumption on the 5p/p 2 momentum resolution and the beamstrahlung 
spectrum. Accuracies are given for an integrated luminosity of 1 ab _1 . 



5p/p 2 


Beamstrahlung 


Fit Result (GeV) 


0. 


none 


1150 ± 10 


3.0 xl(T 5 


none 


1150 ± 12 


4.5 xl(T 5 


none 


1151 ± 12 


4.5 xl(T 5 


Std. 


1143 ± 18 



for the determination of squark masses [7], has been extended also to sleptons for 
the Lhc and a TeV-class LC [8]. It is interesting to consider it here also for its 
implications on the required momentum resolution in the detector. Two values of 
the solenoidal magnetic field B = 4 and 6 T have been tested, corresponding to 
momentum resolutions 5p/p 2 of 4.5 x 10~ 5 and 3.0 x 10 -5 GeV -1 respectively. No 
appreciable difference on the resulting mass accuracy has been observed from these 
two momentum resolution, since the dominant smearing is due to beamstrahlung. 
In fact, at Clic, the main issue is the significant beamstrahlung smearing of the 
luminosity spectrum, and thus of the effective Eb eam value. The corresponding effect 
has been estimated by using both a perfectly well known and constant beam energy 
and that corresponding to the baseline Clic parameters at a nominal y / s=3 TeV. 
Results are summarised in Table 3.1. The smuon mass has been extracted by a x 2 fit 
to the muon energy spectrum by fixing M^o to its nominal value (see Table 3.1. The 
fit has been also repeated leaving both masses free and performing a simultaneous 
two-parameter fit. The results are = (1145±25) GeV and M x <, = (652±22) GeV 
(see Figure 3.1). 



3.2 The Threshold Scan method 

An alternative method to determine the fi mass is an energy scan of the rise of the 
e + e~ — > cross section, close to its kinematical threshold. It has been shown 

that an optimal scan consists of just two energy points, sharing the total integrated 
luminosity in equal fractions and chosen at locations optimising the sensitivities to 
the /2 width and mass respectively [9]. Including the beamstrahlung effect, induces 
a shift of the positions of the maxima in mass sensitivity towards higher nominal y/s 



6 



International Workshop on Linear Colliders 



Table 3: Accuracies on the determinations of the smuon masses for the two bench- 
mark scenarios considered in this study using threshold scans with different experi- 
mental conditions. 



Point 




Beam- 


Pol. 






5M 






strahlung 




(TeV) 


(ab" 1 ) 


(GeV) 


H 


H 


none 


0/0 


3.0-3.5 


1 


± 11 


H 


H 


Std. 


0/0 


3.0-3.5 


1 


± 15 


E 


H 


none 


0/0 


3.8-4.2 


1 


± 29 


E 


H 


Std. 


0/0 


3.8-4.2 


1 


± 36 


E 


H 


none 


80/60 


3.8-4.2 


1 


± 17 


E 


H 


Std. 


80/60 


3.8-4.2 


1 


± 22 



energies. For benchmark point E, the cross section at y/s=3 TeV is too small for an 
accurate measurement. Higher centre-of-mass energy, 4 TeV, and polarised beams 
need to be considered. By properly choosing the beam polarisation, not only the 
pair production cross sections are increased but also their sensitivity to the smuon 
masses. Results are summarised in Table 3. 



4 Scalar Top Analysis 



In the context of the benchmark point E, the process e + e~ — > t\t\, t^ti followed by 
the decay ti.2 — ► tg offers a possibility to determine the scalar top mass based on the 
reconstructed top energy only, following a procedure similar to that adopted for the 
smuon mass. The same relation as Eq. 1 between the top energy spectrum endpoints 
and the masses and yfs energy holds in this case, with M^ and Mg replacing M^ and 
Mo respectively. By determining the endpoints of the top quark energy spectrum 
produced in t — > tg, the masses of the stop and the gluino can be extracted. However, 
here the situation is complicated by the fact that top quarks may also be produced 
in g decays or other t decays, such as t — ► tx°- Figure 4 shows the distribution of 
the energy of top quarks produced in e + e~ — > £2*2 events. Although the expected 
endpoints are diluted, it is still possible to extract M^ from the shape of the top 
energy spectrum. A \ 2 fit) including the contributions from other top sources in 
signal events, indicates that a relative accuracy SMt 2 /Mt 2 = ±7.5% can be achieved 
with 3 ab -1 , without including beamstrahlung effects. 
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Figure 2: Energy spectrum of top quarks in e + e~ — > £2^2 events. The £2 - * tg 
signal is represented by the lower histogram. Left: the different contributions to 
the spectrum. Right: Et op spectrum corresponding to the nominal mass and 3 ab _1 
(points with error bars) compared to the expectations for M t ~ = 1300±60 GeV. The 
endpoint from ti — > tg remains visible when overlayed to the distributions from the 
other decay channels. 

5 Conclusions 

A preliminary study of sfermion reconstruction at the Clic multi-TeV collider has 
shown that smuon masses can be measured with typical accuracies of 1.3-3.0 % for 
C = 1 ab _1 for two benchmark points, using either the muon energy technique, 
when y/s >> 2 x Mx, or a threshold scan. The accuracy of the threshold scans is 
significantly enhanced by the availability of polarised beams. The energy spectrum 
technique has also been extended to the scalar top in the case when the decay ti — * tg 
is kinematically allowed. 

We are grateful to G. Blair, A. De Roeck, J. Ellis, S. Kraml and L. Pape for 
contributions. 
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